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Combining scattering matrix formalism with non-linear cr-model technique we analyze weak local- 
ization effects in arrays of chaotic quantum dots connected via barriers with arbitrary distribution 
of channel transmissions. With the aid of our approach we evaluate magnetoconductance of two 
arbitrarily connected quantum dots as well as of A'' x M arrays of identical quantum dots. 
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I. INTRODUCTION 

Quantum interference of electrons is fundamen- 
tally important for electron transport in disordered 
conductorsiiSi^. Quantum coherent effects are mostly 
pronounced at low temperatures in which case certain 
interaction mechanisms are effectively "frozen out" and, 
hence, cannot anymore restrict the ability of electrons 
to interfere. At the same time, there exists at least one 
mechanism, electron-electron interactions, which remains 
important down to lowest temperatures and may destroy 
quantum interference of electrons down to T = 0. It 
is, therefore, highly desirable to formulate a general the- 
oretical formalism which would allow to describe elec- 
tron interference effects in the presence of disorder and 
electron-electron interactions at any temperature, includ- 
ing the hmit T ^ 0. 

In a series of papers^ we offered such an approach 
which extends Chakravarty-Schmid description'^ of weak 
localization (WL) and generalizes Feynman- Vernon path 
integral influence functional technique to fermionic sys- 
tems with disorder and interactions. With the aid of this 
technique it turned out to be possible to quantitatively 
explain low temperature saturation of WL correction to 
conductance 5G^^{T) commonly observed in diffusive 
metallic wires^. It was demonstrated'* that this satura- 
tion effect is caused by electron-electron interactions. 

It is worth pointing out that low temperature satura- 
tion of WL correction and of the electron decoherence 
time (extracted from 5G'*^'"{T) or by other means) 
has been repeatedly observed not only in metallic wires 
but also in virtually any type of disordered conductors 
ranging from individual quantum dots^ to very strongly 
disordered 3d structures and granular metals'^. It is quite 
likely that in all these systems we are dealing with the 
same fundamental effect of electron-electron interactions. 
In order to support (or discard) this conjecture it is nec- 
essary to develop a unified theoretical description which 
would cover essentially all types of disordered conductors. 
Although the approachi^ is formally an exact procedure 
treating electron dynamics in the presence of disorder 
and interactions, in some cases, e.g., for quantum dots 
and granular metals, it can be rather difficult to directly 
evaluate 6G'^^{T) within this technique. 

One of the problems in those cases is that the de- 



scription in terms of quasiclassical electron trajectories 
may become insufficient, and electron scattering on dis- 
order should be treated on more general footing. An- 
other (though purely technical) point is averaging over 
disorder. In our approach** disorder averaging is (can be) 
postponed until the last stage of the calculation which is 
convenient in certain physical situations. In other cases 
- like ones studied below - it might be, in contrast, more 
appropriate to perform disorder averaging already in the 
beginning of the whole analysis. In addition, it is desir- 
able to deal with the model which would embrace various 
types of conductors with well defined properties both in 
the long and short wavelength limits. 

In this paper we make a first step towards this uni- 
fied theory. Namely, we will describe a disordered con- 
ductor by means of an array of (metallic) quantum 
dots connected via junctions (scatterers) with an arbi- 
trary distribution of transmissions of their conducting 
channels. This model will allow to easily crossover be- 
tween the limits of granular metals and those with point- 
like impurities and to treat spatially restricted and spa- 
tially extended conductors within the same theoretical 
framework. Electron scattering on each such scatterer 
will be treated within the most general scattering ma- 
trix formalismJiS adopted to include electron-electron 
interaction cffectsiSiiiiiSiiiii. Averaging over disorder 
will be performed within the non-linear cr— model tech- 
nique in Keldysh formulation as first proposed by Hor- 
bach and Schoni^ for non-interacting electrons. This 
method has certain advantages over the imaginary time 
approach*^ since it allows to treat both equilibrium and 
non-equilibrium problems and also enables one to include 
Coulomb interaction between electrons in a straightfor- 
ward manner*''. In this paper we will analyze WL correc- 
tions to conductance merely for non-interacting electrons 
and will include interaction effects only phenomenologi- 
cally by introducing an effective electron dephasing time 
Tip as a parameter of our theory. Systematic analysis of 
the effect of electron-electron interactions on weak local- 
ization within this formalism will be developed elsewhere. 

The structure of the paper is as follows. In Sec. II 
we define the basic model of a Id array of quantum dots 
and outline the key features of our formalism. In Sec. 
Ill we will introduce convenient parameterization of the 
non-linear cr-model which will then be used in Sec. IV 
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to derive WL correction to the system conductance for 
the model in question. This WL correction will be eval- 
uated for various structures in Sec. V. Sec. VI contains 
direct generalization of our analysis and results to the 
case of 2d arrays of quantum dots and is followed by a 
brief summary in Sec. VII. Some technical details of our 
calculation are presented in Appendix. 



II. THE MODEL AND FORMALISM 

Let us consider a Id array of connected in series chaotic 
quantum dots (Fig. QJ. Each quantum dot is character- 
ized by its own mean level spacing (5„. Adjacent quan- 
tum dots are connected via barriers which can scatter 
electrons. Each such scatterer is described by a set of 
transmissions of its conducting channels Tjf'^ (here k la- 
bels the channels and n labels the scatterers). We will 
ignore spin-orbit scattering and, for the sake of definite- 
ness and simplicity, we will first focus our attention on Id 
arrays only. Generalization of our analysis to other situ- 
ations can be performed in a straightforward manner, as 
it will be demonstrated in Sec. VI of the paper. 

An effective action S[Q] of an array depicted in Fig. 
depends on the fluctuating 4x4 matrix fields^* Qn{ti, ^2) 
defined for each of the dots {n — 1, ...,N — 1). Each of 
these fields is a function of two times ti and t2 and obeys 
the normalization condition 



1. 



(1) 



The action of an array can be represented as a sum of 
two terms 



iS[Q] = iSd 



iSt 



(2) 



The first term, iSd[Q], describes the contribution of bulk 
parts of the dots. This term reads 



iSdlQ] 



a, 



(3) 





FIG. 1: Id array of A'^— 1 quantum dots coupled by barriers. 
Each quantum dot is characterized by mean level spacing S„ . 
Each barrier is characterized by a set of transmissions of its 

(n) 

conducting channels T^. . 



A similar expression was also considered within the imag- 
inary time techniqueiSiSi. 

Note that here the magnetic field H is included only 
in the term Q describing the quantum dots while it is 
ignored in the term |(SJ). Usually this approximation re- 
mains applicable at not too low magnetic fields. We will 
return to this point in Sec. VI. 

An equilibrium saddle point configuration A{ti — 12) 
of the matrix field Q{ti,t2) depends only on the time 
difference and has the form 
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(6) 



-g'^iE) -1 ; 



where g^{E) = 2[1 - 2/f(S)] = 2 tanh(i;/2r). This 
choice of the saddle point corresponds to the following 
structure of the 4x4 matrix Green function G: 
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Here we defined the time inversion operator T: 
Tf{t) = f{tf~t), 



(7) 



(8) 



where t/ will be specified later. Note that the function 
G in Eq. Q, defined for a given disorder configuration, 
should be contrasted from the Green function 



Here H is an external magnetic filed, an = 
bn{e^ /fi'^c?)vFdn min{^e, 6„ is a geometry dependent 
numerical prefactoriSiiS, dn is the size of n— th dot, le is 
the elastic mean free path in the dot, and A is 4 x 4 
matrix: 
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The second term in Eq. (|2Jl, iSt[Q], describes electron 
transfer between quantum dots. It has the forr^ 
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defined for a given realization of the matrix field Q. In 
Eq. lO we also introduced the electron elastic mean free 
time Tp. 



III. GAUSSIAN APPROXIMATION 

In order to evaluate the WL correction to conductance 
we will account for quadratic (Gaussian) fiuctuations of 
the matrix field Qn- This approximation is always suffi- 
cient provided the conductance of the whole sample ex- 
ceeds e^//i, in certain situations somewhat softer applica- 
bility conditions can be formulated. Expanding in powers 
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of such fluctuations we introduce the following parame- 
terization 

= A + i[Wn,k] + WnkW,, ~ \{W^,,k} + 0{W^){10) 

It follows from the normalization condition (QJ that only 
8 out of 16 matrix elements of W are independent pa- 
rameters. This observation provides certain freedom to 
choose an explicit form of of this matrix. A convenient 
parameterization to be used below is 



I Mi„ hxn 

U2„ h2n 

ain + bin Wi„ 

V a2n + b2n V2n 



(11) 



With this choice the quadratic part of the action takes 
the form 



(12) 



(2) 

where iS^j^ [a, b] does not depend on H and describes dif- 

(2) 

fuson modes, while iSuJ [u, v] is sensitive to the magnetic 
field and is responsible for the Cooperons. The diffuson 

(2) 

part of the action iS^^f [a, b] was already analyzed in Ref. 
[ij and will be omitted here. Below we will focus our 
attention on the Cooperon contribution which reads 



iS^J [u, v] 



n—1 
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at 



n—1 

y ^Tr 

n=l 



dt 



[v2n,Vin] - 16a„ff^ V1V2 



{Uln - Ul.n-l){u2n " W2,n-l) 



+ (win - Wl,n-l)(w2n - W2,«-l) 



(13) 



where g„ = 2^j,T^"^ — 27r?i/e^i?„ is the dimensionless 
conductance of n— th barrier. With the aid of the action 
((T^ we can derive the pair correlators of the fields ui^2 
and V12' 

{Uln{tl,t2)u2m{t',t")) = {vin{t' , t" )v2rn{h , t2)) 

= ^^(il -t2+t'- nCnn^it" - h), (14) 

where we defined a discrete version of the Cooperon 
Cnm{t) obeying the equation 
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[{gn + gn+l)Cn 



- gn+lCn + l,m] = 5nm5{t). (15) 



This equation should be supplemented by the bound- 
ary condition Cnm{t) = which applies whenever one 



of the indices n or m belongs to the lead electrode. Here 
THn — l/16a„7f^ is the electron dephasing time due to 
the magnetic field. In Eq. p5|l we also introduced an 
additional electron decoherence time in n— th quantum 
dot Tipn which can remain finite in the presence of in- 
teractions. In this paper we are not aiming to further 
specify the interaction mechanisms and only account for 
them phenomenologically by keeping the parameter 
in the equation for the Cooperon. 



IV. WL CORRECTIONS 

Let us now derive an expression for WL correction to 
the conductance in terms of the fluctuating flelds u and v. 
In what follows we will explicitly account for the discrete 
nature of our model and specify the WL correction for a 
single barrier in-between two adjacent quantum dots in 
the array. 

We start, however, from the bulk limit, in which case 
the Kubo formula for the conductivity tensor ctq/j reads 

(Jai3{r,r') I dt' {t-t') 

J — OG 

X (jpit',r')j^it,r) - j„it,r)jpit',r')). (16) 

Following the standard procedur o^i^i'^ , approximating the 
Fermi function as —dfF{E)/dE « 5{E) (which effec- 
tively implies taking the low temperature limit) and using 
a phenomenological description of interactions as medi- 
ated by external (classical) fluctuating fields, from Eq. 
(|16|l one can derive the WL correction in the form: 



(5(7 



WL, 



r,r') 



dt' I dt 



X (G«(t,ri;i",r^)G-4(t',ri;t,r2)),.^^ 



max cross 



(17) 



which implies summation over all maximally crossed 
diagrams^-^-'^, as indicated in the subscript. At the same 
time, averaging over fiuctuations of Q within Gaussian 
approximation is equivalent to summing over all ladder 
diagrams. Since we are not going to go beyond the above 
approximation, we need to convert maximally crossed di- 
agrams in Eq. IjlTfl into the ladder ones. Technically this 
conversion can be accomplished by an effective time re- 
versal procedure for the advanced Green function which 
can be illustrated as follows. 

Consider, e. g., the second order correction to in 
the disorder potential Udis{x) 

(5(2)G"4(i',rJ;t,r2) = -■>' dT2 ^ "^^^ / ^^^^d^ici 

X C^it', r[;Ti,Xi)Udis{xi)G^{Ti,Xi;T2,X2) 
xUdisix2)G^{T2,X2;t,r2). (18) 
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Making use of the property G^(Xi, X2) = G^*{X2,Xi), 
we get 



X G^*(i, rz; T2, a;2)[/dis(a;2)G^*(r2, 0:2; n, a;i) 
xC/di,(a;i)G«*(ri,a;i;t',r;). (19) 

Setting tf ~ t + t' , we rewrite this expression as follows 

Jtf-t Jtf-t 
X J d^X2d^xi G^*{tf ~t',r2;T2,X2) 

X Udm{x2)G^*{T2,X2]Ti,Xi) 

xUdisixi)G'^*iTi,Xi;tf^t, r[). 



(20) 



Close inspection of the right hand side of Eq. H2()|l allows 
to establish the following relation 

5(2)G^(t',rl;i,r2) -r<5(2)G«*(t',r2;t,rl)r, (21) 

which turns out to hold in all orders of the perturbation 
theory in Udis- As before, the time inversion operator T 
is defined in Eq. (jSJ with tf =t + t'. 

As a result, the expression for Sa^J" takes the form: 



'5</(r,r') 



dt' / dt 



)r\ —r%—r' 



X (G^(i,ri;t",ri)TG^*(i',r2;t,ri)T) 



dis, ladder 



(22) 



Rewriting Eq. (|22ll in terms of the matrix elements of 
the Green function Q, we obtain 



dt' / dt" 



X (G33(i,ri;t",r^)G44(t',r2;t,r;)),., 1^,^^^ (23) 

Our next step amounts to expressing WL correction 
via the Green function Gq . For that purpose we will 
use the following rule of averaging 

(G33(i,n;t",ri)G44(i',r2;t,ri)),i^ 

= (G33;Q(i, ri;t'\ r'^)Gii,Q{t\ r2; t, r^))^ 

- (G34;Q(t,ri;t,ri)G43;Q(t',r2;t",r^))Q. (24) 

One can check that within our Gaussian approximation 
in u and v the first term in the right hand side of Eq. 
(|24|l does not give any contribution. Hence, we find 



dt' j dt" 



X (G34;Q(t,ri;t,ri)G43;Q(t',r2;i",ri))Q. (25) 



Let us now turn to our model of Fig. 1 in which case 
the voltage drops occur only across barriers. In this case 
Eq. p5(l . which only applies to bulk metals, should be 
generalized accordingly. Consider the conductance of an 
individual barrier determined by the following Kubo for- 
mula 



G 



dt'{t-t'){I{t',x')I{t,x) 



I{t,x)I{t',x')). 



(26) 



Here I{t, x) is the operator of the total current flowing 
in the lead (or dot) and a; is a longitudinal coordinate 
chosen to be in a close vicinity of the barrier. Due to 
the current conservation the conductance G should not 
explicitly depend on x and x'. Comparing Eqs. 126() 
and pt)|l . and making use of Eq. (|25|l and the relation 
I{t, x) — J d^zjxit, X, z), where jx is the current density 
in the a;— direction and z is the vector in the transversal 
direction, we conclude that WL correction to the conduc- 
tance of a barrier between the left and right dots should 
read 



SGYi 



dt' / dt" / d'^zd^z' 



Aum? 

X Xi X2)xi—X2—x{^ x'-^ x'^x'^—x'^—X' 

X (G34;Q(t, Xi , Z] t, x'l , 2;')G43;Q (*', ^2, Z; t" , x'2, z')) ^ . 

(27) 

In what follows we will assume that both coordinates 
X and x' are on the left side from and very close to the 
corresponding barrier. Let us express the Green function 
in the vicinity of the barrier in the form 



GQ{t,x,z;t',x',z') 



Y^{e^P'^^^-^P'^-'g++(t,t',x,x') 



+ e^P-^+^P-^'g+-{t,t',x,x') 



g-+{t,t',x,x')}^^izW^{z'), (28) 



where ^n{z) are the transverse quantization modes 
which define conducting channels, p„ is projection of 
the Fermi momentum perpendicular to the surface of the 
barrier, and the semiclassical Green function slowly 
varies in space. Eq. H27I) then becomes 



47rTO2 



dt' / dt" 



E 



E 

mnkl a/375— ±1 



X e 



{Qli,,,(t,t,x,x')gilS'^t''^^^^'))^ 

ia.pnX-i—i(3p.n-^x'-^ -\-i'ypkX2—iSpi x'^ 



(29) 



Next we require 5G^^ to be independent on x and 
x', i.e. in Eq. H29|) we omit those terms, which contain 
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quickly oscillating functions of these coordinates. This 
requirement implies that ap^ + ^Pk = and (3pm + ^Pi — 
0. These constraints in turn yield 7 — —a, S — — /3, 
k = n and I = m. Thus, we get 



E E 

mn q/3=±1 



dt' / dt"a(3p 

nPn 



Q 



(30) 



Let us choose the basis in which transmission and 
reflection matrices t and f are diagonal. In this ba- 
sis the semiclassical Green function is diagonal as well, 
Qmn oc Gnn^nrm and Eq. H3U|I takcs the form 
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dt' / df 



x(^tnn;34(i>i)^Z,;„;43(i'.^") 



' L,nn;34 



++ 

L,nn;43 



(t'X) 



^L,rm;34(^' ^)^L,nn;43(* ' * ) 



Q' 



(31) 



What remains is to express WL correction in terms of 
the field Q only. This goal is achieved by establishing an 
explicit relation between the Green function Q and the 
field Q. A derivation of this relation is presented in Ap- 
pendix A. Here we only display the final result expressed 
via the fluctuating fields vi and V2- We obtain 



6G 



WL 
LR 
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TT 



E 



dt' / dt" 



X {Tr,[viL{t,t)v2R{t',t")+ViR{t,t)v2L{t',t")] 

+ T^iviLit, t) - ViRit, t)] [v2Lit', t") - V2R{t', <")])• (32) 

Note that the contribution linear in T„, which contains 
the product of the fluctuating fields on two different sides 
of the barrier, vanishes identically provided fluctuations 
on one side tend to zero, e.g. if the barrier is directly 
attached to a large metallic lead. In contrast, the contri- 
bution cx r,j in Eq. H32I) survives even in this case. 
Finally, applying the contraction rule (|I4|I we get 



5GY^ 



9 

eg 

"47r2 



dt{(3[5RCLR{t)+5LCRL[t)\ 



{I- P)[5RCRR{t)+5LCLL{t)\]. (33) 



Here 5l^r is the mean level spacing in the left/right quan- 
tum dot, g — '2J2k'^k is the dimensionless conductance 
of the barrier and j3 = 71- (I — T^)/ T^. is the cor- 
responding Fano factor. 

Likewise, the WL correction to the n— th barrier con- 
ductance in Id array of A'^ — 1 quantum dots with mean 
level spacings Sn connected by N barriers with dimen- 
sionless conductances f/„ and Fano factors /3„ reads 



SG. 



WL 



47^2 



dt{Pn[SnGn-l,n{t) 




FIG. 2: Single quantum dot connected to the leads via two 
barriers. 



+ Sn-lGn^n-lit)] + (1 " f^n) [SnGnnit) 

+ ,5„_iC„_i,„_i(t)]}. (34) 

So far we discussed the local properties, namely WL 
corrections to the conductivity tensor, Sa^^{r,r'), and 
to the conductance of a single barrier, SG^r. Our main 
goal is, however, to evaluate the WL correction to the 
conductance of the whole system. For bulk metals one 
finds that at large scales the WL correction H17() is lo- 
cal, Sa^0 {r,r') oc S{r ~ r'). In general though, there 
can exist other, non-local, contributions to the conduc- 
tivity tensor—. Without going into details here, we only 
point out that, even if these non-local terms are present, 
one can still apply the standard Ohm's law arguments 
in order to obtain the conductance of the whole sample. 
Specifically, in the case of Id arrays one finds (cfi^): 



5G 



WL 



1 



1 



higher order terms. (35) 



Eqs. (123), ESI and will be used to evaluate WL 
corrections for different configurations of quantum dots 
considered below. 



V. EXAMPLES 
A. Single quantum dot 

We start from the simplest case of a single quantum 
dot depicted in Fig. 2. In this case the solution of Eq. 
(|15|l reads 



Cii [t) = exp 



t t 



(36) 



where tjj = Air / [gi + g2)Sd is the dwell time, and 5d is 
the mean level spacing in the quantum dot. All other 
components of the Cooperon are equal to zero. From 
Eq. we get 
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FIG. 3: Most general system with two quantum dots 

According to Eq. the total WL correction becomes 
WL 5i.9l(l- A)+5?52(l-/32) 



5G 



47r2 (51 + 32)2 (I/td + l/r^ + \/th) ' 



(38) 



bmce I/th oc i?^ the magnetoconductance has the 
Lorentzian shaped. In the limit H — and in the ab- 
sence of interactions {r^p ^00) Eq. H38() reduces to24 

.WL _ giglil - Pi) + gh2{l - P2) 



6G' 



(39) 



71" {gi+g2f 

B. Two quantum dots 

Next we consider the most general setup composed of 
two quantum dots with the corresponding conductances 
and Fano factors defined as in Fig. |3| The Cooperon 
is represented as a 2 x 2 matrix which zero frequency 
component satisfies the following equation 



511 + 312 + gy +71 - 9v 

-gy 521 + 522 + 5y + 72 
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" I 471/62 



where 
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(41) 



Defining A = igii+gi2+gy + 7i){92i+g22+gy + l2)-gy, 
we get 



gy/^i (.911 +512 +5y + 7i)/^2 



With the aid of Eq. (|33|) we derive WL corrections for 
all five barriers in our setup: 
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FIG. 4: The magnetoconductance of two dots of Fig. |^for 

di,d2 » le, di/d2 = 5, Qij = go, fiij = 0, /3y = 0, r^i = 
Tip2 — 00. Here Hi = l/Ay/aiToi is the field at which weak 
localization is effectively suppressed in the first dot. For gy = 
the magnetoconductance is given by superposition of two 
Lorentzians with different widths (decoupled dots), while for 
large gy only one Lorentzian survives corresponding to the 
contribution of a one "composite dot" . 
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(42) 



WL correction to the conductance of the whole struc- 
ture 8G^^ is obtained from the general expression for 
the conductance determined by Ohm's law: 

G = [GiiGi2(G2i + G22) + G2iG22(Gii -I- G12) 
4- Gy{G\2 + G22)(Gii -I- G21)] 
/[(Gii-HGi2)(G2i-f G22) 
+ Gj,(Gii + Gi2 + G2i + G22)]. (43) 



Substituting Gy — > Gy -|- SG^J^ into this formula and 



FIG. 5: Two quantum dots in series. 



expanding the result to the first order in SG^^, we get 



— 1,2 ^ 



(44) 



Combining Eqs. (|42|I - H44() we arrive at the final result 
for the WL correction to the conductance of the whole 
structure. This general result is rather cumbersome. It 
is illustrated in Fig. 21 for a particular choice of the sys- 
tem parameters. Below we will specifically consider two 
important limits. 

First we analyze the system of two quantum dots con- 
nected in series, as shown in Fig. [S] i.e. in the general 
structure of Fig. |2|we set G12 — G21 — 0, Gn = Gi, 



G^ 



G2, G22 = G3, = (iy = jii and fi^i = /^s 



We also assume H = and r^p ~ 00. WL corrections to 
the barrier conductances then take the form 



SGI 



WL 



WL 



SG. 



gi(g2 + .g3)(i- A) 

5152 + 5253 + 5i53 ' 

g2(gi+g3)(i-/92) + 2ffi 

gi52 + aigj, + gig3 
53(gi +g2)(i-/33) 



5152 + 5253 + 5i53 
while Eq. H43|) reduces to 

G1G2G3 



G 



G\Gi + G1G3 + G2G3 
WL correction for the whole system then reads 



(45) 



(46) 



SG 



WL 



5i.9i5i(.92 +53)(l-/?i) 
(5152 + 5253 + 5153)^ 

5i5253(gi +g3)(l -fe) 
1" (51.92 + 5253 + 5153)^ 

5i5l53(.9i +.92)(1 -fe) 
I" (5152 + 5253 + 5153)^ 
2e^ 5i5i5i 

TT (5152 + 5253 + 5153)^ ' 



(47) 



In the limit of open quantum dots, i.e. Pi,2,3 — 0, we 
reproduce the result of Ref. It is easy to see that 
provided the conductance of one of the barriers strongly 
exceeds two others, Eq. (|T7|) reduces to Eq. If 
all three barriers are tunnel junctions, /3i,2,3 1, the 
first three contributions in Eq. (|47|l vanish, and only 
the last contribution - independent of the Fano factors - 
survives in this limit. If, on top of that, one of the tunnel 



FIG. 6; The system of two connected quantum dots only one 
of which is in turn connected to the leads. 

junctions, e.g. the central one, is less transparent than 
two others, 52 <C 51,53, the result acquires a particularly 
simple (non-Lorentzian) form 



2e2 



9l 



TT (51 +71) (.93+72)' 



(48) 



with 71,2 defined in Eq. (gTJ. Note that SG^^ cx 55, 
this result is dominated by the second order tunneling 
processes across the second barrier. 

Our second example is the system depicted in Fig. [HI 
which corresponds to the following choice of parameters 
in Fig. El Gn = Gi, G12 = G2, G21 = G22 = 0, /3n = /3i 
and /3i2 — lii- In addition, we assume that electrons are 
subject to dephasing only in the second quantum dot, 
i.e. r^i = CX) while r^2 is finite. This setup allows one 
to analyze the so-called dephasing by voltage probe3^^*2&. 
We obtain 



Gn 



47r 



5y + 47r/(52r^2 



S\ (51 + 52).9y + 47r(5i + 52 + gy)IS2T^2 ' 



In the limit T^p2 00 this result reduces to 

47r 1 



Gn = 



5i 51 + 52 



(49) 



(50) 



and we again arrive at Eq. (|39|l . i.e. the second quantum 
dot attached to the first one does not affect the expres- 
sion for WL correction. In the opposite limit of short 
decoherence times, t^2 0, we find 



Gn = 



47r 



1 



Si 91+92+ gy 

and arrive at the WL correction2& 

5i.9i(l-/3i)+5?52(l-/32) 



SG^L^ 



^ i9l+92)Hl + TD/Tf) 



where 



9y 



51 + 52 TD 



(51) 



(52) 



(53) 



is the electron decoherence rate induced in the first quan- 
tum dot due to coupling to the second one acting as an 
effective voltage probe. 
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C. ID array of identical quantum dots 

Let us now turn to Id arrays of quantum dots depicted 
in Fig. ^ For simplicity, we will assume that our array 
consists of A'' — 1 identical quantum dots with the same 
level spacing 5„ = Sd and of N identical barriers with the 
same dimensionless conductance gn = g and the same 
Fano factor (3n = f3. We will also assume that the quan- 
tum dots have the same shape and size so that thu = th 
and Tipn = . For this system the Cooperon can also be 



found exactly. The result reads 



2 

c„„i(cj) = — ■ 

9=1 



■nan ■ Trom 

sm sm -j^ 



1 

Th 



1-cos ^ 



(54) 



Here to = 2iT/g5d and th = \/lQaH^. 
The WL correction then takes the form 



5G 



WL 



^1 /3cos^ + l-/3 



9=1 



(55) 



The sum over q can be handled exactly and yields 



6G 



WL 



ttN^ 



N- 



1 — u 



2N 



(3(1 + u^) + 2{1 - (3)u 



(7V-l)/3 



where 




(56) 



(57) 



oo our ri 



esult 



In the tunneling limit P — 1 and for r^p 
defined in Eqs. H56() - (|57l) becomes similar - though not 
exactly identical - to the corresponding resulliSi. 

If is long enough, namely l/r^ -Exh , where i?Th = 
tt'^ /2N^Tu is the Thouless energy of the whole array, in 
Eqs. H55|l - H56(l it is sufficient to set — oo. In this 
case the magnetic field H significantly suppresses WL 
correction provided 1/th ^ i?Th or, equivalently, if 



1 / irgSd 



H Hat, II at — ,, 

~ ' 8N\I a 

In the opposite limit l/r^p > i?Th we find 



(58) 



SG 



WL 



e 



1 + ID + la 

TH Ty, 



1 



(59) 



In particular, in the diffusive limit th^t^ ^ tu we get 
-,WL _ / Dtht^ 



irNd y Th +T, 



(60) 




FIG. 7: Magnetoconductance of a Id array of A'^ — 1 identical 
open (/3 = 0) quantum dots in the absence of interactions 
(r^ (x). The field Hn is defined in Eq. 1581 . 



where we introduced the diffusion coefficient 



D = (f/2TD. 



(61) 



Eq. H60|l coincides with the standard result for quasi- 
Id diffusive metallic wire. Note, however, that the val- 
ues of Th within our model may differ from those for a 
metallic wire. The ratio of the former to the latter is 
''''h/'''h°^ ~ T-fl/Tz), where ta ~ d/vp is the flight time 
through the quantum dot. Since typically ta < td we 
conclude that for the same value of D the magnetic field 
dephases electrons stronger in the case of an array of 
quantum dots. 

For a single quantum dot {N = 2) Eq. (|56|) reduces to 



SG 



WL 



e2(l-/3) 



1 



An 


('^ 


TD 






TH 





(62) 



in agreement with Eq. H38|l . 

For two identical quantum dots in series we obtain 



5G 



WL 



e 



2~p 



1 



2to 
Th 



2to 



1 



2td 
3th 



2to 
3r„ 



, (63) 



i.e. the magnetoconductance is just the sum of two 
Lorentzians in this case. 

Finally, in the absence of any interactions {t^p = oo) 
and at = we obtain 



6G 



WL 



e 



I 
N 



1 

iV2 



(64) 



In the limit N ^ oo this result reduces to the standard 
one for a long quasi-ld diffusive wireSi while for any finite 
iV we reproduce the results for tunnel barriersSi (/3 — > 1) 
and open quantum dotsSS (/3 ^ 0) . 
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FIG. 8: 2d array of identical quantum dots. Here the number 
of barriers in the a;-direction is chosen to be A'' = 4, and the 
number of quantum dots in the y— direction is M = 2. The 
barriers are characterized by the dimensionless conductances 
Qx and Qy and the Fano factor (3. 



VI. GENERALIZATION TO 2D ARRAYS 

Until now our analysis was only focused on structures 
with several quantum dots and Id arrays. Generaliza- 
tion to the case of 2d and 3d systems is straightforward. 
Below we analyze an important case of 2d arrays. 

Consider an array consisting of — 1 x M quantum 
dots. For simplicity, here we will only deal with the case 
of identical quantum dots, see Fig. 7. The WL correction 
to the conductance of this array reads 



N M 



6G 



WL 



WL 

nrn ' 



(65) 



where, similarly to Eq. H33() 

WL _ ^^9x^d 



5G. 



47r2 

+ C'n.ri-l 



(ii|/3 [C„_i,„;mm(0 

{t)] + {I - l3)[Cnn-mAi) 

m (66) 



defines the WL correction for the barrier with "coordi- 
nates" n,m. In order to find the Cooperons Cnn';mm'{t) 
one needs to solve the equation 

5 1 1 \ „ 5^ 



47 [(^^^ 



2.g«) 



^ dt TH 

^ Cnn' \nini' 9xCn—l^n'-.mra' 9xGri+l^n' \rava' 

9yCn,n' ~l\mm' QyCn.n' + l\mm''\ — f^nn' t^mm' ^(^) . (67) 

which is directly analogous to Eq. (|15|l . In the zero 
frequency limit the solution of this equation with appro- 
priate boundary conditions reads 

Cnn';mm'(w > 0) — YTat ^ ! ^ ^ 



MN 



■ N 



M 



■iTqy(m' — 1) 



M 



7b + i + ^(l-cosl^) + ^(l-cosl^ 



N 



COS 



irqyVa 
M 



-\- COS 



■n-i?y(m— 1) 



M 



S0,q, + 1 



(68) 



Combining Eqs. (|^ . and we obtain 



SG 



WL 



a 9x 



„ 7V-1 M- 

27r27V2 2^ 

9x = l <?H=0 



E 



/3cos^ 



(69) 



This result is valid provided an external magnetic field H 
dephases electrons predominantly inside quantum dots. 
This is the case provided the field is not too low'^*'-'^^, H > 
H* , where H* = TrcgSra/ed? . At lower fields, H < H* , 
one can apply the standard theory^ developed for homo- 
geneous metals, in which case l/rf^'^^ — AeDH/c, where 
D is now defined in Eq. H61|) . Substituting the value 
H* into this expression we get l/r|f ~ tr/tI). Compar- 
ing this energy scale with and i?Th we immediately 
arrive at the condition 



max{l/r^,£;Th} > 5^<5^Tfl, 



(70) 



for which Eq. H69|) is applicable at all values of H. 

Turning to concrete examples we first consider the sim- 
plest case with N = 2,M = 2, which is a symmetric 
version of the system of Fig. 13 Eq. (|69|l then yields 



6G 



WL 



1 



e'^gxSd 



1-/3 



87r2 J_ + J_ 



ga:Sd 

2tt 



27T 



For H = and oo this expression reduces to 



g^WL ^ _ 



e2(l-/3) 



An 



1 



9x + gy 



Next we consider an extended isotropic {gx 
array and stick to the diffusive regime 



(71) 



(72) 



2d 



Em 



In this case we find 



< 



1 



SG 



WL 



e2 M 
'2^ iV 



In- 



TH 



TD 



td{th 



2.773- 7r/3 



.(73) 



M 



The leading term in this equation matches with the stan- 
dard WL correction for a 2d diffusive metallic film in the 
parallel magnetic field2. 

Let us briefiy discuss an effect of anisotropy. In the 
limit of small gy <C gx/N"^ the system reduces to a set of 
M essentially independent Id arrays and, hence, SG^'" = 
MSGf/^, where 6Gf/^ is defined in Eqs. |HSEni- In the 
opposite limit of large gy ^ gxM"^ electron diffusion in 
the direction perpendicular to the current becomes fast, 
and one can treat the system as a Id array of — 1 
composite quantum dots each of them consisting of M 
original dots. In this limit we get 6G^^ = SG^/^ . 
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Finally, let us note that our Eq. ifBU]! also allows to 
reproduce recent results^ for WL correction to the con- 
ductivity of bulk granular metals. In order to handle this 
limit, in Eq. Ht)9|l one should formally set M, N ^ oo 
(which yields SG^^ oc M/N and allows to define the 
conductivity) and then put P = \ and gx = Sy- 



VII. SUMMARY 

In this paper we have developed a theoretical approach 
based on a combination of the scattering matrix formal- 
ism with the non-linear cr-model technique. This ap- 
proach allows to analyze weak localization effects for an 
arbitrary system of quantum dots connected via barri- 
ers with arbitrary distribution of channel transmissions. 
This general model can be used to describe virtually any 
type of disordered conductors. Employing our approach 
we have evaluated WL corrections to the system con- 
ductance in a number of important physical situations, 
e.g., for the case of two quantum dots connected to each 
other and to external leads in an arbitrary way (Sec. V 
B), as well as for Id (Sec. V C) and 2d (Sec. VI) arrays 
of identical quantum dots. In a number of specific lim- 
its our general results reduce to those derived earlier by 
means of other approaches. 

The results obtained here remain valid either in the 
absence of interactions or provided the interaction ef- 
fects on weak localization are taken into account within a 
phenomenological scheme which amounts to introducing 
electron decoherence time Ti^ as an additional parameter. 
The method proposed here also serves as a good start- 
ing point for a more general and systematic analysis of 
electron-electron interaction effects. This analysis will be 
worked out in our forthcoming publications. 

This work is part of the EU Framework Programme 
NMP4-CT-2003-505457 ULTRA-ID "Experimental and 
theoretical investigation of electron transport in ultra- 
narrow 1-dimensional nanostructures" . 



Defining the diagonal matrix v = VnSnrm and making 
use of the normalization condition we can write the 
solution in the form 



e«"^/2-<=(l - «Qi) + e--"-/2-=(l + aQL) 
-ia5af3V~^d{x — x') -f E^l^ 

fi-i.'/2re(i -f aQL) + e-«"'"'/2-e(l ^ uQl) .(A2) 



Here i?^'^ is an arbitrary operator. Requiring not to 
grow exponentially far from the barrier we arrive at the 
following constraints: 



1 + Ql 
1-Ql 



Rl Rl 
1-Ql 
1 + Ql 



0, 
= 0. 



(A3) 



Similarly, for the right dot we obtain 
1~Qb. \ f RV R 







1 

Rt 
Rr" Rr 



R 



^R 



1 



Rr Rr 
-Qr 
1-Qr 



0, 
= 0. 



(A4) 



Note that the elastic mean free time Tg drops out of Eqs. 
(|A3IA4|I . thus indicating a very general nature of these 
constraints. The Green functions on the left and right 
barrier sides are related to each other by the 5*— matrix 



f i' 
i f' 



of this barrier. This relation has the form 



(A5) 



(A6) 



APPENDIX A: RELATION BETWEEN GREEN 
FUNCTION g AND Q. 



M 



r f't'-' 



(A7) 



Here we will closely follow the method proposed by 
Nazaro\«Sli. Let us select one of the barriers in our ar- 
ray and denote (coordinate-independent) Q— fields in the 
left (right) dot with respect to this barrier as Ql {Qr)- 
Provided Ql.r are slow functions of time, in the barrier 
vicinity one can neglect the term id/dt. In addition, one 
can linearize the electron spectrum in the vicinity of the 
Fermi energy and replace V^/2m ±ivmd/dx, where 
Vn — Pn/rn is the electron velocity in a given channel. 
As a result, for the left dot one gets 



^ d I - 



Qlt = Six - x')6„ 



,6af3. (Al) 



being the transfer matrix which satisfies M<TzM^ = <Jz- 

Eqs. (|X3|) . l|X4|l and (|X6l) for R°^jf^ can be resolved 
making use of the fact that in the barrier vicinity, i.e. for 
I a; I, ^ WnTe, the Green function takes the form 



^L,R. — 



Rtji 

Rl%. 



R 
R', 



■L,R 







■L.R 

sign(a;i - 0:2) - Ql,r _ 1 
sign(xi - X2) - Ql,r ' ' 



The operators R^^j^ turn out to be diagonal in the chan- 
nel indices in the basis for which the matrices t and f 
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are diagonal as well. Defining the channel transmission 
values r„ — |i„P, we get 



(All) 



'•L,nm 



i + ^({QL,Qfl}-2) 



T„ {qr + [Qr,Ql]~QlQrQl) _ tI{1-Ql) 

8 2 I (A9) 

t„{1+Ql) r„ (QR-[Qn,Q^]-QLQflQi. j ' 
2 8 



a/3 
R^nm 



1 + ^{{Ql,Qr}-2) 



t„[Ql+IQr. 



JrQlQr) r;(l+QR) 
2 



2 8 



(AlO) 



In order to find WL correction H3U|) it is sufficient to 
determine the Green function ljA9p only in the left dot: 



^L.nm 



R 



2v„ 



L,7im 



2Vn 



( 2Aui 0\ 


5^Aui 
\ -Ai;25^ 2Ai)2 ) 

/ \ 

-2Au2 

Aui5^' -2Aui 

y -g-^Aua y 



where Aui,2 = ui,2fl — ui,2L, Awi,2 — fi,2fl — 'yi,2L and 



-2i 

2U2L 



-2uiL \ 



-viLg^ ~ g^uiL 2viL 



\ v2Lg + g 



K 



-2v2L -2i J 



Here R'^'^m has been expanded to the first order in ui,2 
and vi^2- From Eqs. IjASp and (jAll|) we find 

C++ ^— _-t'iL-T„Ai;i 



(7 + + _ V2L+TnAv2 Q 



Substituting these expressions into Eq. H31|) . after some 
transformations we arrive at Eq. (EH). 



L,nn;34 ~ 




L.nn; 34 


VlL 

~rn 

Vn 




V2L 


L,nn ;43 


Vn ' 


-+ - 


V2L 


L,nn;43 


r-a • 

Vn 
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